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ABSTRACT 

We present radio transient search algorithms, results, and statistics from the ongoing Arecibo Pulsar 
ALFA (PALFA) survey of the Galactic plane. We have discovered seven objects through a search for 
isolated dispersed pulses. All of these objects are Galactic and have measured periods between 0.4 
and 4.7 s. One of the new discoveries has a duty cycle of 0.01%, smaller than that of any other radio 
pulsar. We discuss the impact of selection effects on the detectability and classification of intermittent 
sources, and compare the efficiencies of periodicity and single-pulse searches for various pulsar classes. 
For some cases we find that the apparent intermittency is likely to be caused by off-axis detection or 
a short time window that selects only a few bright pulses and favors detection with our single-pulse 
algorithm. In other cases, the intermittency appears to be intrinsic to the source. No transients 
were found with dispersion measures large enough to require that they originate from sources outside 
our Galaxy. Accounting for the on-axis gain of the ALFA system, as well as the low gain but large 
solid-angle coverage of far-out sidelobes, we use the results of the survey so far to place limits on the 
amplitudes and event rates of transients of arbitary origin. 

Subject headings: 



1. INTRODUCTION 

Radio pulsars show a wide variety of modulations 
of their pulse amplitudes, including bursts and nulls, 
that affect their detectability in surveys. Phenom- 
ena seen in some pulsars include short-period nulling, 
in which a pulsar is not detected for severa l pulse 
perio ds, only to reappear with full strength ([Backed 
Il970f ): eclipses, in which a companion star or its wind 
or magnetosphere absorbs or disperses the pulsar sig- 
nal (Fruchter et al. 1988l lStappers et al. 1996. Lvn e et al.l 
[r993.Kaspi et al. ,2004) : long-term nulling or intermit- 
tent beha vior, in which a pu lsar is quiescent for days 
or weeks (jKramer e t al."20 Q6|); and rota ting radio tran- 
sients (RRATs. iMcLaughhn et al.l[2006h . pulsar-like ob- 
jects from which only occasional radio bursts are de- 
tected. This paper describes analysis of a large-scale 
survey using the Arecibo telescope that is sensitive to 
both periodic and aperiodic signals. 



RRATs were firs t discovered in archive P arkes Multi- 
beam survey data ([McLaughlin et al] 120061 ). Eleven ob- 
jects, with periods ranging from 0.7 to 7 s and pulse 
widths of 2 — 30 ms, were found u sing a single pulse 
search algorithm ([McLaughlinl l2007f ) . The longer peri- 
ods of RRATs compared with the general pulsar pop- 
ulation suggest similarities with the X-ray populations 
of X-ray dim isolated neutron stars (XDINSs) and mag- 
netars. RRAT J1819-1458 has b een detected at X-ray 
energies ([McLaughlin et all 120071 ) with properties that 
are similar to those of XDINSs and high magnetic field 
radio pulsars. 

A different type of pulse modulation is observed in 
the case of pulsars emitting giant pulses. Such pulses 
are tens to thousands of times brighter and an order 
of magnitude or m ore narrower than the average pulse 
fsee iKnightj |2006| for an overview). Giant pulses from 
the Crab pulsa r have substructure on timescales of 2 ns 
([Hankins et all 120031 ) . and PSR B1937+21 emits giant 
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pulses as narrow as 16 ns (|Popov et al.ll2004f) . Giant 
micro-pulses from the Vela pulsar have widths ~ 50/xs 
(jJohnston et al.l [200 11). and the slowly rotating pulsars 
PSR B1112+50, PSR B0031-07, and PSR J1752+2359 
occasionally emit bright pulses which are 1 — 10 ms wide, 
5 — 3 times narrower than the average pulse ijKnightl 
I2006D . The detection of giant pulses is a potentially pow- 
erful method for finding extragalactic pulsars too distant 
for their normal emission to be detect able by periodicity 
searches (jMcLaughlin fc Cor dcs 20031. 

A variety of energetic phenomena other than pul- 
sar emission can give rise to fast transients potentially 
detectable in radio pulsar surveys. Within the Solar 
System, transient radio events may be generated by 
energetic particles impacting the Earth's atmosphere, 
solar flares, and decameter radio flares originating in 
Jupiter's atmosphere. Analogously to the latter, extra- 
solar planets with strong magnetic fields ar e expected to 
be detectable in the 10 - 1000 MHz range (|Farrell et all 
119991 iLazio et al.l |2004 iZarka et all I2001D . Magnetic 
activity on the surfaces of brown dwarfs and particle 
acceleration in the magnetic fields of flare stars are 
also known r adio flare progenitors (iBerger "eTall [2nnil 
Bergeij [20021 iGarcia-Sanchez et al.l l2003l iJackson et al.l 
1989t). Gamma-ray bursts are predicted to have de- 
tectable radio emission at ~ 100 MHz dU sov fc Katd 
I2OOOI ISagiv fc WaxmaiJl2002[ ) , and radio flares have been 
observed from som e X-ray binaries (jWaltman et al.|[r995l 
iFender et al.lll997f l. Among the most energetic and ex- 
otic events in the Universe, supernovae, merging neu- 
tron stars and coalescing black holes may produce wide- 
band radio bursts detecta ble at extragalactic distances 
([Hansen fc Lvutikovl[200l . 

In this paper we describe an ongoing survey for pulsars 
and transient radio sources with the Arecibo telescope. 
The survey addresses outstanding questions about the 
nature and emission mechanisms of intermittent radio 
sources. In § [5] we present the PALFA survey parame- 
ters, and in §[3] we describe the single pulse search and 
radio frequency interference excision algorithms which 
are part of the survey data processing pipeline. Sec- 
tion m contrasts PALFA detection statistics on known 
pulsars and new discoveries, and § [5] examines selec- 
tion effects influencing the detection and classification 
of transient sources. In § [HI we apply an intermittency 
measure method for comparing the efficiency of period- 
icity and single pulse pulsar searches. In § [71 we discuss 
the properties of individual intermittent objects discov- 
ered by PALFA. In § [51 and § [9l we apply constraints 
derived from the survey sensitivity and results to the de- 
tectability of various energetic phenomena expected to 
emit radio bursts. Finally, in § llOi we present our main 
conclusions. 

2. PALFA SURVEY OBSERVATIONS 

2.1. Survey Parameters 

The PALFA survey started in 2004, shortly after the 
installation of the seven-beam ALFA receiver on the 
Arecibo telescope. The survey searches for pulsars and 
transients in the inner and outer Galactic plane regions 
accessible to Arecibo (see below). 

The ALFA receiver is well-suited for survey observa- 
tions, allowing simultaneous data collection from seven 
fields, each 3.5' (FWHM) across. Taking into account 



the hexagonal arrangement of the beams on the sky and 
the near sidelobes, th e combined power- pattern is ap- 
proximately 24' X 26' ([Cordes et al.ll2006[ ). We observe 
a 100 MHz passband centered on 1440 MHz in each of 
the seven telescope beams. Wid eband Arecibo Pulsar 
Processors (WAPPs; iDowd et"aLi,2000) are used to syn- 
thesize 256-channel filterbanks spanning these bands at 
intervals of 64 /xs. During observations, full-resolution 
data are recorded to disk and, in parallel, decimated 
down to a time resolution of 1 ms and searched for peri- 
odic signals and single pulses by a quick-look processing 
pipeline running in real time at the Arecibo Observatory 
([Cordes et al.l[2006[ ) . This approach allows for immediate 
discovery of relatively bright pulsars with periods longer 
than a few milliseconds. Searching full-resolution data al- 
lows detection of millisecond pulsars and narrower single 
pulses and is done offline at participating PALFA insti- 
tutions as the processing is much more computationally 
intensive. 

Table [1] lists various ALFA system and survey param- 
eters, including the sky area corresponding to processed 
and inspected data reported on in the pres ent paper and 
the t otal sky area observed to date (see iGordes et al.l 
I2OO6I for a detailed explanation of other parameters). 
Standard observation times are 268 s for inner Galaxy 
pointings (30° < I < 78° ,\b\ < 5°) and 134 s for outer 
Galaxy pointings (162° < I < 214°, |6| < 5°). Some early 
observations had a duration of 134 s and 67 s for inner 
and outer Galaxy pointings, respectively. The quoted 
system temperature of 30 K is measured looking out of 
the Galactic plane. The initial threshold of 5<t for the 
single pulse search is used when selecting events from 
dedispersed time series based on their signal-to-noise ra- 
tios only, before any filtering is applied. While there 
are a significant number of events due to random noise 
above this threshold, identification of an event as a gen- 
uine pulse takes into account not only its signal-to-noise 
ratio but also the fact that it is detected at a contiguous 
range of trial dispersion measures, which is in general 
not true of spurious events. Thus weak pulses can be 
correctly identified, while they would be excluded if the 
threshold was set according to Gaussian noise statistics. 

2.2. Survey Sensitivity 

Here we compute the maximum distance at which 
sources of a given luminosity are be detectable by the 
PALFA survey, I?max- We then compare the sensitiv- 
ity of the PALFA survey to previous w ork done with 
the Parkes Multib eam system (Manches ter et al.l |2001L 
iMcLaughlin et al.ll2006f l. 

The rms noise in a radio transient search, where the 
effective observation time is equal to the pulse width, is 



^iVpoi A/ W' 



(1) 



where Sgys is the system-equivalent flux density, A^poi = 2 
is the number of polarization channels summed and A/ 
is the bandwidth. A pulse's observed width W may be 
broadened compared to its intrinsic width Wi by sev- 
eral effects. After dedispersing the raw data and obtain- 
ing a dedispersed time series, there is residual dispersive 
broadening due to the finite width of a frequency channel 
and the error of the trial dispersion measure (DM) used 
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PALFA SURVEY PARAMETERS. 




Parameter 


Value 


Center frequency (GHz) 


1.440 


Total bandwidth (MHz) 


100 


Channel Bandwidth (MHz) 


0.39 


Sampling time (fis) 


64 


Nominal system temperature (K) 


30 


Gain (K/Jy) 




Center pixel 


10.4 


Ring pixels 


8.2 


Beam width (') 




Main beam (1 pixel) 


3.5 


Main beam plus near sidelobes (1 pixel) 


18 


Main beam plus near sidelobes (7 pixels) 


24 X 26 


Inner Galaxy 




Observation time per pointing (s) 


268, 134 


Observed (deg^) 


156 


Processed^ (deg^) 


99 


Anticenter 




Observation time per pointing (s) 


134, 67 


Observed (deg'^) 


119 


Processed^ (deg^) 


87 


Detection threshold (cr) 




Single pulse search (initial threshold) 


5 


Fourier Transform periodicity search 


7.5 
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^ Full-resolution data processed with the Cornell pulsar 
search code. 

compared to the actual pulsar DM. Scatter broadening 
is not correctable and will have a contribution that de- 
pends on observing frequency and varies with direction 
on the sky. In general, 
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where AtDM,ch — 8.3//sDMA/ch//'^ is the dispersive 
broadening across a frequency channel of width A/ch 
(MHz) for an observing frequency / (GHz), AtDM,crr is 
the dispersive broadening due to the difference between 
the trial and actual DM of the source, and Atsc oc f^^ 
is the scattering broadening. Broadening conserves pulse 
area, so that the intrinsic and observed peak flux densi- 
ties are related through Sp \ Wi 



SpW. If 5p,„i„ 



is the detection threshold, the minimum detectable in- 
trinsic peak flux density is 
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For a one steradian pulsar radio beam, a source of intrin- 
sic peak luminosity Lp.i can be detected out to a maxi- 
mum distance of 
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(4) 

For a steadily emitting pulsar with period-averaged lu- 
minosity L and duty cycle /dc, we have Lp « L/ f^c- 

The amount of pulse broadening depends on sys- 
tem parameters as well as dispersion and scatter- 
ing, which vary with direction on the sky so that 
W = W{Wul,bJ,Af,N,i,,Ssys,At). We use the 
NE2001 model of G alactic ionized electron density 
(jCordes fc Lazidl2002t ) to calculate representative results 
for i^max in the direction I — 35°, 6 = 0°, a region of 
overlap between PALFA and the Parkes Multibeam sur- 
vey. Fig. [1] shows I?max VS. Lp^i detection curves using a 



threshold m = 6 for both surveys. For lower luminosities, 
sources are not visible to large enough distances for scat- 
tering to affect detectability and the inverse square law 

1/2 

dominates the detection curve so that Dmax oc Lp ; . For 
larger distances and smaller intrinsic pulse widths, scat- 
tering and (residual) dispersion smearing make pulses 
increasingly harder to detect and Umax increases more 
slowly with ip,i. 
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Fig. 1. — Maximum distance at which transients with vari- 
ous peak luminosities and intrinsic pulse widths of (top to bot- 
tom for each set of curves) 50, 5, and 0.5 ms can be detected by 
the PALFA and Parkes Multibeam surveys. The linear portion 
of the curves corresponds to a luminosity-limited regime. In that 
regime, Dmax is about twice as large for PALFA than for the Parkes 
Multibeam survey because of the larger sensitivity of the Arecibo 
telescope. Breaks in the curves correspond to a transition from 
luminosity-limited to scattering-limited detection for PALFA and 
to a dispersion-limited regime for Parkes. The difference is due 
to the smaller channel width of PALFA (0.4 MHz) compared to 
Parkes (3 MHz). 



3. SINGLE PULSE SEARCH METHODS 

This section describes processing methods employed by 
the Cornell pulsar search pipeline^ at the Cornell Cen- 
ter for Advanced Computing and the Swinburne Uni- 
versity of Technology. The PRESTO search code ^ is 
run independently at West Virginia University, Univer- 
sity of British Columbia, McGill University and the Uni- 
versity of Texas at Brownsville. PRESTO uses a similar 
matched filtering algorithm to the one described below, 
but a different RFI excision scheme and trial DM list. 

We dedisperse raw data with 1272 trial DMs in the 
range 0—1000 pc cm~^. In order to find individual pulses 
from intermittent sources we operate on dedispersed time 
series with two time domain algorithms: matched filter- 
ing, which has been the standard in single pulse searching 
so far, and a friends-of-friends algorithm. After single 
pulse candidates have been identified we use a stack- 
ing method in the time-frequency plane to verify that 

^ http : //arecibo .to. Cornell . edu/PALFA/ 
^ http: //www . cv .nrao . edu/~sransoni/presto 
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the pulses are dispersed and the sweep observed across 
frequencies follows the dispersion relation, as expected 
for non-terrestrial sources. Certain types of terrestrial 
signals, e.g. swept-frequency radars, have pulses whose 
appearance in the time-frequency plane may approxi- 
mate dispersion by ionized gas, but on closer inspection 
such pulses generally deviate significantly from the cold 
plasma dispersion relation. 

3.1. Matched Filtering 

Matched filtering detection of broadened pulses relies 
on the convolution of a pulse template with the dedis- 
persed time series. Ideally, a pulse template would con- 
sist of one or several superimposed Gaussian templates, 
reflecting the diversity of pulsar pulse profiles, some of 
which exhibit multiple peaks. In addition, multipath 
propagation due to scattering adds an exponential tail to 
the pulse profile. We approximate true matched filtering 
by smoothing the time series by adding up to 2" neigh- 
boring samples, where n = — 7, and selecting events 
above a threshold after each smoothing iteration. The 
smoothing is done in pairs of samples at each n-stage, so 
that the resulting template is a bo xcar of length 2" sam- 
ples (jCordes fc McLaughlinI [20031 ) . The samphng time 
used in PALFA observations is 64 fj,s and therefore our 
matched filtering search is most sensitive to pulse widths 
of 64 IIS to 8.2 ms. This search strategy is not optimal 
for single pulses from heavily scattered pulsars because 
the templates are symmetric, while the scattered pulse 
shape with an exponential tail is not, and significantly 
scattered pulses can be wider than ~ 10 ms. 

3.2. Finding Time-Domain Clusters 

Two issues make a complement to matched filtering 
necessary. The pulse templates described above have 
discrete widths of 2" samples by algorithm design and 
there is decreased sensitivity to pulses with widths that 
are significantly different. In addition, the matched fil- 
ter search output may be dominated by bright, wide 
pulses from radio frequency interference (RFI). In that 
case, a single RFI burst is detected as an overwhelm- 
ing number of individual events instead of a single event. 
The cluster algorithm is similar to the friends-of-friends 
search algorithm used t o find galaxies in optical images 
(jHuchra fc Gellei1ll982l ). and complements matched fil- 
tering by not restricting the width of expected pulse de- 
tections. The dedispersed time series is processed se- 
quentially. Each event above a threshold which is found 
is designated as the first of a cluster. A cluster of events 
is augmented by broadening it to include any adjacent 
samples above a threshold. Gaps of Jigap samples are al- 
lowed within a cluster, and PALFA data is processed with 
'^gap = 2. The brightest sample of a cluster is recorded 
as the event amplitude, and the total number of samples 
in the cluster as its width. This approach is less sensitive 
to weak, narrow pulses but results in significantly fewer 
spurious events due to RFI. 

A limiting factor for the largest pulse width detectable 
by both the matched filtering and friends-of-friends 
search algorithms is the fact that data are analysed in 
blocks much shorter than the complete time series length. 
The mean and standard deviation of a block are used for 
thresholding events found within the block. This ap- 
proach minimizes the effect of baseline variations with 
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Fig. 2.— Single pulse search output for PSR J0627+16. The 
bottom panel shows events with S/N > 5 vs. time and DM with 
larger circles denoting brighter bursts. Panels on top from left 
to right show histograms of the number of events vs. S/N and 
DM, and a scatter plot of event DM vs. S/N. A fit to the nar- 
row peak of event DM vs. S/N indicates a pulse width ~ 2 ms 
llCordes fc McLaughlin|[200a) . 

time scales much longer than typical pulsar pulse widths 
of a few to a few tens of ms. The disadvantage is that 
only pulses with W ^ Tbiock can be detected. In pro- 
cessing PALFA data, we use blocks of length Tbiock = 
4096 X 64 ^s = 0.26 s. Accord ing to the NE2q01 mo del 
of ionized gas in the Galaxy (jCordes fc Laziol [2002( 1. a 
scattering broadening time of that magnitude in the 
inner-Galaxy part of the PALFA survey corresponds to 
DM > 1000 pc cm""^ and a maximum search distance 
well outside of the Milky Way for most directions we 
survey. However, lines of sight intersecting HII regions 
can result in large dispersion measures and scattering 
times and are therefore selected against. Any intergalac- 
tic scattering that broadens pulses beyond about 0.1 s 
will also cause events to be selected against. 

3.3. Example Results 

Fig. [5] shows standard single pulse search output for 
the discovery observation of PSR J0627-hl6 (§ 8.1). The 
main panel shows signal-to-noise ratio (S/N) vs. DM 
and time for events with S/N > 5. The panels on top 
illustrate event statistics for the observation: from left 
to right, number of events vs. S/N and DM, and DM 
vs. S/N. The pulsar detection is manifested as a peak 
in the histogram of number of events vs. DM. There 
is a corresponding peak in the DM vs. S/N plot, since 
the several detected pulses are significantly brighter than 
background noise in the dedispersed time series. 

Fig. [3] shows events in DM-time space for another 
PALFA single pulse discovery, PSR J1928-hl5 (§ 8.5). 
In this case, three closely spaced bursts were found at 
t ~ 100 s, DM - 250 pc cm-3. The pulsar is detected 
at a range of DMs, with the signal to noise ratio increas- 
ing as trial DM values approach the actual pulsar DM 
and decreasing as trial DMs further on recede from the 
pulsar DM. In contrast, events due to interference from 
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Fig. 3. — In the upper panel, events with S/N > 4 from the 
discovery PALFA observation of PSR J1928+15 are displayed in 
the DM-time plane. Events aligned vertically and spanning all 
trial DM values at t ^ 46,71,82,90,118 s are due to terrestrial 
interference. Three closely spaced bursts were found at t ~ 100 s, 
DM ~ 250 pc cm"''. The lower panel shows a magnification of 
the area around the bursts in the dedispersed time series (top) and 
the DM-time plane (bottom). The intervals between successive 
bursts are 0.403 s, establishing the pulsar period. Larger circles 
denote brighter bursts, and the brightest burst has S/N ~ 60 in 
the dedispersed time series. The scaling of circle size with signal- 
to-noise ratio is slightly different in the two plots. 

terrestrial signals a,t t ^ 46, 71, 82, 90, 118 s span the en- 
tire range of trial DMs and their signal to noise ratios do 
not show a significant variation with DM. 

If an excess of candidate pulses is identified in the 
dedispersed time series for a particular trial DM, we use 
the expected dispersion sweep across the frequency band 
in order to test if the pulses are from non-terrestrial ori- 
gin. For a pulse with S/N > 10 in the time series, we 
extract a chunk of raw data centered on its time of ar- 



FlG. 4. — Dispersion sweep in the time-frequency plane of the 
brightest single pulse detected from (a) PSR J1928-I-15 and (b) 
PSR J1946-I-24; (c) Stacked dynamic spectrum of the five brightest 
pulses detected from PSR J0627-f 16; (d) Stacked dynamic spec- 
trum of the two brightest pulses detected from PSR J1909-f06. 
The pulses are dispersed, such that the higher frequencies arrive 
earlier. 

rival and look for a sweep across frequency that follows 
the dispersion relation (Fig. [4] a, b). 

3.4. Time- Frequency Plane Stacking 

In the case of multiple weaker pulses detected in the 
same beam, we extract a raw data chunk centered on the 
arrival time of each pulse, stack the chunks and look for 
a dispersion sweep in the resulting cumulative dynamic 
spectrum (Fig.[4]c, d). 

Stacking can induce spurious dispersion sweeps in the 
time-frequency plane even when the data consist of only 
Gaussian-distributed noise. Single pulse events identi- 
fied in dedispersed time series by definition correspond 
to sums over a given dispersion path that are above av- 
erage. When these are used to select chunks for sum- 
ming, the statistical fluctuations will build up to show a 
dispersed pulse in the time-frequency plane that follows 
the cold-plasma dispersion law perfectly. We simulated 
this effect by generating a time-frequency plane of Gaus- 
sian noise, dedispersing it, and selecting events above 
threshold from the resulting time series. We extracted 
chunks from the fake time-frequency plane centered on 
each event and stacked them, thus reproducing the pro- 
cedure used on real data. 

Selecting events with signal-to-noise above a threshold 
TUd in the dedispersed time series yields a biased set of 
noise-only samples in the time-frequency plane. On the 
other hand, when stacking dynamic spectra, points in the 
stacked dynamic spectrum must have a minimum signal- 
to-noise TTic 2 for a dispersion sweep to be discernible 
by eye. The rms noise in the dedisper sed t ime series after 
summing Nch channels is a factor of ^/N^h larger than the 
rms noise in the time-frequency plane. The rms noise in 
the stacked dynamic spectrum after stacking chunks 
is a factor of \/W^ larger than the rms noise in the time- 
frequency plane. We calculate the average deviation from 
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the mean in the time-frequency plane imphed by the two 
thresholds and rric, equate them, and obtain 

Nc > {mjmdf iVch. (5) 

For the PALFA survey and our processing pipeline, 
iVch = 256 and rrid = 5. Therefore if Nc > 40 data 
chunks are added there will be a spurious, dispersed pulse 
in the resulting stacked dynamic spectrum even if the 
chunks contain only noise. Since the noise decorrelates 
over 1 — 2 samples, the pulse will have a width on the 
order of 64—128 /iS. To avoid the induced spurious-event 
effect, the width and S/N of stacked dynamic spectra 
should be compared with those that would result from 
noise only. All of our detections that use stacking satisfy 
the criteria W >> 64 fis and Nc << 40. Typically we 
sum chunks centered on the brightest 2 — 5 pulses for 
a single pulse candidate, which is safely below the Nc 
hmit. In the case of PSR J1909-h06 (Fig. Hi) only 2 
chunks were added and the dispersed swath in the time- 
frequency plane is 1 ms wide. 

3.5. RFI Excision 

The RFI environment at the Arecibo telescope and 
the 7-beam configuration of the ALFA receiver present 
both challenges and opportunities for RFI mitigation 
to facilitate searching for single pulses. PALFA survey 
data are dedispersed with trial dispersion measures of 
— 1000 pc cm~^. RFI pulse intensity typically peaks 
at DM = pc cm^^, and incidental low-intensity pulses 
whose S/N peaks at DM = pc cm~'^ tend to be smeared 
below the detection threshold for dedispersed time series 
with trial DM > 50— 100 pc cm~^. A more complex sig- 
nature is observed for Federal Aviation Administration 
(San Juan airport) radar pulses, which are unfortunately 
common in Arecibo observations. The radar rotation pe- 
riod is Pr = 12 s, and each pulse has an envelope that 
is ~ 1 s wide and consists of sub-pulses with variable 
period on the order of 2 — 3 ms. Depending on telescope 
orientation, radar pulses may be detected in all, some, 
or none of the ALFA beams due to reflections off the 
telescope structure. Radar pulses are up to two orders 
of magnitude brighter than pulsar pulses and, without 
mitigation, can completely dominate single pulse search 
results. In addition, the modulation of the radar signal is 
manifested as detections with DM 7^ pc cm~'^ so that 
unlike other non-radar RFI, their S/N vs. DM signature 
cannot be used for excision. 

We exploit the known radar characteristics as well as 
the pattern of pulse detection in the 7 ALFA beams to 
excise both radar and non-radar RFI. The first part of 
our excision algorithm targets radar pulses. After a list of 
single pulse events is generated for all trial DMs, we bin 
the events for trial DM = — 3 pc cm~'^ in time (by 0.1 s) 
and record the number of events in each time bin. Then 
we treat the histogram as a time series and perform a dis- 
crete Fourier transform. A peak near the radar rotation 
frequency indicates a significant number of radar pulses 
in the data. From the Fourier components we extract 
the phase and find the arrival time of the earliest radar 
pulse, to- Since the envelope width is ~ 1 s, events within 
0.5 s of that location are excised, and the procedure is 
repeated for events near t = to -\- NPr, where N is an in- 
teger. The more radar pulses are present, the better the 
performance of this technique because the radar peak in 



the DFT is more prominent and the pulses' arrival times 
are determined more precisely. However, if only a couple 
of strong radar pulses are present within the typical 268 s 
PALFA integration time, they may be bright enough to 
dominate event statistics, yet the DFT method does not 
excise non-periodic incidental RFI. Consequently, after 
applying the DFT-based method we use an additional 
RFI filter that handles aperiodic cases. 

The second filter uses the number and proximity of 
beams in which an event is detected in order to deter- 
mine if it is due to RFI. Again events for trial DM — 
— 3 pc cm~^ from each beam are binned in time. Af- 
ter detecting peaks indicating an excess of events for the 
respective time bins, the algorithm cross-checks between 
results for all seven beams, and each event falling within 
a histogram peak receives a penalty grade based on how 
many beams' histograms exhibit a peak and how close 
to each other they are on the sky. Most pulsars de- 
tected blindly via a single pulse search appear in one 
beam or two adjacent beams, and very bright pulsars 
may be detected in three or four adjacent beams. We set 
the excision penalty threshold just below the value cor- 
responding to the latter configurations and excise events 
accordingly. This method complements the DFT clean- 
ing scheme and effectively excises sparse radar blasts as 
well as non-periodic RFI detected in multiple beams. 
The application of the two excision algorithms makes a 
marked difference in the final single pulse search output 
for pointings contaminated with RFI (Fig. . The fig- 
ure shows some false positives, for example a clump of 
pulsar pulses in beam 4 around t — 80 s are excised due 
to low- level RFI at low DMs occurring in non- adjacent 
beams in the same 0.1 s time bins as the pulses. Lower- 
ing the threshold according to which an excess of events 
is defined and time bins are marked for excision reduces 
false positives but also diminishes the effectiveness of RFI 
excision. While tens of pulses are detected within 134 s 
from the bright pulsar B2020-I-28 (Fig.O, for sources dis- 
covered via a single pulse search the number is an order 
of magnitude smaller (see § 7). Therefore the chance of 
RFI occurring simultaneously with pulses of intermittent 
sources and causing them to be excised is much lower. 

4. PALFA PULSAR DETECTION STATISTICS 

The PALFA survey has made a total of 354 blind de- 
tections of 172 pulsars up to late 2008. Following quick- 
look processing of the data with degraded time and fre- 
quency resolutions at the time that data are acquired 
(jCordes et al.ll2006[ ). data are shipped to external sites 
for processing through the Cornell pulsar search code 
and PRESTO. Table [H shows a breakdown of Cornell 
detections, including pulsars originally discovered by the 
quick-look pipeline, by the periodicity (FFT) and sin- 
gle pulse (SP) search algorithms for known pulsars and 
new discoveries, and Table[2]lists parameters for PALFA 
single pulse discoveries that are discussed in more detail 
in § [71 Columns in Table [3] are (1) Pulsar name; (2)- 
(3): equatorial coordinates; (4) period; (5) pulse width 
(FWHM); (6) dispersion measure; (7) distance; (8) peak 
flux density of the brightest detected pulse; (9) total ob- 
servation time; (10) total number of pulses detected; and 
(11) implied average pulse rate. The last column (12) 
points to the paper section that discusses the object. 

More than half of the detected known pulsars and only 
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Fig. 5. — Single pulse search output for a blind detection of 
pulsar B2020+28 before (top) and after (bottom) excising radar 
and incidental RFI. Each row shows results from one ALFA beam. 
From left to right, panels show: events with SjN > 5 vs. DM 
channel and time, number of events vs. DM channel, and event 
S/N vs DM channel. The pulsar detection in beams 3 and 4 is 
evident after RFI events are excised. 



TABLE 2 

Pulsar detection statistics by algorithm. 



Pulsars 


EFT only 


SP only 


EFT and SP 


Total 


Known 


48 (38%) 


5 (4%) 


73 (58%) 


126 


New 


28 (61%) 


6 (13%) 


12 (26%) 


46 



a quarter of the newly discovered objects were seen by 
both the periodicity and single pulse search algorithms. 
There is also a significant difference in the fraction of 
pulsars detected only via a single pulse search. In total, 
74% of the new pulsars as opposed to 42% of the known 
pulsars were detected either only by FFT search or only 



by single pulse search. The much higher percentage of 
new pulsars to be detected only by one algorithm means 
that most of the PALFA discoveries are either periodic 
emitters too weak to show up in single pulse searches or 
intermittent objects whose emission is too heavily mod- 
ulated be detected by a periodicity search. Considering 
that bright pulsars are more likely to be seen by both 
the FFT and SP search, and they are also more likely to 
have already been found by previous, less sensitive sur- 
veys, this indicates that PALFA is probing deeper and 
finding pulsars that are farther away or less bright than 
the known objects in the same region of the Galactic 
plane. 

Single pulse searches have not been routinely car- 
ried o ut on survey data in the past, with the excep- 
tion of iPhinnev fc Tavloii (fTOTOf) . iNicel (fl999V who dis- 
covered PSR J19 1 8-I-08 via a single pulse search, and 
iMcLaughlin et al.l (|2006f ). who found 11 RRATs. There- 
fore, the proportion of single-pulse-only detections is 
higher for the PALFA discoveries than the known pul- 
sars. The five known pulsars detected only by single 
pulse search are known steady emitters which were seen 
away from the beam center and therefore with reduced 
sensitivity. 

5. SELECTION EFFECTS AND INTERMITTENCY 

The sample of PALFA single pulse discoveries (§[7]) in- 
cludes one object which has not been detected again de- 
spite reobservations (J1928-I-15), two sources detectable 
in re-observations through their time-averaged emission 
(PSR J0628-I-09 and PSR J1909+06), and four long- 
period objects, one of which has the smallest known duty 
cycle for any radio pulsar (PSR J0627-I-16). 

Among the 11 Parkes RRATs, 10 were successfully 
redetected in subsequent Parkes observations, and one 
(J1839— 01) has not been redetec ted despite multiple at- 
tempts (|McLaughhn et al.ll2006f) . PSR J0848-43 and 
PSR J1754— 30 can sometimes be detected through their 
time-avera ged emission with s ensitive, low-frequency ob- 
servations ( McLaughlinI |200"7[) . Six of the sources have 
periods greater than 2 s. 

In this section we discuss selection effects which may 
account for some of the observational characteristics of 
intermittent sources. Characteristics of this population 
which remain unexplained by observational biases may 
indicate underlying intrinsic differences between these 
objects and conventional radio pulsars. They may be- 
long to a genuinely intermittent and physically different 
class of neutron stars. 

5.1. Multiplicative Effects 

We can define a general model for observed signal in- 
tensity in terms of the time and frequency-dependent flux 
density Si{t, v) as 

/ (t, I.) = Gt (t, v) Giss (t, v) Si {t, iy)+n {t, i')+RFl {t, v) , 

(6) 

where Gt is the telescope gain, Gjss is the variation fac- 
tor due to interstellar scintillation, nit, v) is radiometer 
noise, and RFI{t, v) is the contribution from terrestrial 
radio frequency interference. 

Variations in telescope gain across the beams of the 
multi-beam receiver or within the power pattern of a 
single beam impact the detectability of a pulsar as ei- 
ther a single pulse or periodic source. As shown by 
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TABLE 3 

Parameters of PALFA single pulse discoveries 



Pulsar 


RA 


DEC 


P 


W 


DM 


D b 




Ttot 


A^tot 


Rate 


Comment 




(hh:mm:ss) 


(dd:mm) 


(s) 


(ms) 


(pc cm^'^) 


(kpc) 


(mJy) 


(s) 




(h-i) 






J0627+16 


06:27:13(7) 


16:12(2) 


2.180 


0.3 


113 


3.2 


150 


7454 


48 


23 


§ 


7.1 


J0628+09 


06:28:33(7) 


09:09(2) 


1.241 


10 


88 


2.5 


85 


1072 


42 


141 


§ 


7.2 


J1854+03 


18:54:09(7) 


03:04(2) 


4.559 


50 


216 


5.5 


14 


388 


9 


84 


§ 


7.3 


J1909+06 


19:09:24(7) 


06:40(2) 


0.741 


1.5 


35 


2.2 


82 


536 


10 


67 


§ 


7.4 


J1919+17 


19:19:47(7) 


17:44(2) 


2.081 


100 


148 


5.3 


12 


393 


35 


320 


§ 


7.5 


J1928+15 


19:28:20(7) 


15:13(2) 


0.403 


5 


242 


7.4 


180 


2880 


3 


4 


§ 


TE 


J1946+24 


19:46:00(7) 


23:58(2) 


4.729 


4 


96 


4.3 


101 


268 


4 


54 


§ 


7.7 



Position uncertainties eorrespond to the angular radius (out to 50% of boresight power) of an individual ALFA 
beam in tlie discovery observation.'' Estimate from tlie NE2001 model of Galactic eletron density (Cordes & Lazi^ 
nOO^ ." Peak flux density at 1.4 GHz defined as 5 = {S/N)Ssys/ ^/N~[AfW, where S/N is the signal to noise ratio 
of the brightest detected pulse. 



the reobservations of PSR J0627+16, PSR J1909+06, 
PSR J0848-43 and PSR J1754-30, classifying a source 
as intermittent may be due to insufficient sensitivity to 
detect its periodic emission during the discovery obser- 
vation. The effectiveness of single pulse vs. FFT-based 
periodic search depends on the pulse amplitude distri- 
bution, and the observed pulse amplitude distribution in 
turn depends on the pulsar location with respect to the 
telescope beam power pattern. Fig. [S] shows the prob- 
ability density function (PDF) of event amplitude for 
PSR B1933-I-16, a strong pulsar detected in two beams 
of the same PALFA pointing. One beam of the ALFA 
receiver is pointing directly at the pulsar and the low- 
intensity tail of the distribution is clearly visible since 
all pulses are above the detection threshold. The pul- 
sar is also detected in the first sidelobe of an adjacent 
beam, in which case only the brightest pulses are above 
the detection threshold and only the high-intensity tail 
of the distribution is visible. Thus an off-axis detection 
of a steadily emitting pulsar may misrepresent it as an 
intermittent source. 

A different set of effects which affect detectability in 
single pulse vs. periodic search is due to intervening 
ionized gas. Turbulence in the ionized interstellar gas 
and the motion of the pulsar with respect to the gas in- 
troduce phase modulations in the pulsar emission which 
cause the observed intensity to change over time and 
frequency. In general, the scintillation gain factor has 
contributions from both diffractive (DISS) and refractive 
(RISS) interstellar scint illation, such that Giss{t,'^) = 
GDiss{t,v)Gmss{t,v) (|Rickettl[l990l ). The detectabil- 
ity of nearby pulsars with low DMs can be significantly 
affected by DISS, which can have time scales compara- 
ble to typical survey pointing times (minutes) and fre- 
quency scales A/diss << /• In terms of the observing 
frequency / and the pulsar distance D, the DISS time 
scale AiDiss oc The scintillation bandwidth 

^/diss E^nd scattering timescale Tg a re related through 
27r A/diss Ts ~ 1- An empirical fit from lBhat et al.l ()2004l ) 
for / in GHz and Tg in ms gives an estimate for Tg based 
on DM: 

logTg = -6.46+0.154 log DM+1.07(log DM)^-3.86 log/, 

but there is significant scatter about this relationship. 
When the scintillation bandwidth and time scale are 
small, averaging over the observation time and band- 
width can quench DISS. 
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Fig. 6.— Observed single pulse PDFs for PSR B1933+16 in two 
beams of a PALFA pointing: on source (right) and 6 arc minutes 
away (left, shaded). All pulses are detected individually in the on- 
source beam, while only the high intensity tail of the distribution 
is detected via a single pulse search in the off-source beam. A 
sidelobe detection of a canonical pulsar will pick out the brightest 
pulses and may misrepresent the pulsar as an intermittent source. 

For refractive scintillation, A/riss oc / and A^diss oc 
jo.57^,1.6^ Unlike for DISS, the time scale of RISS in- 
creases with distance and therefore it can become very 
long for distant pulsars with relatively high DMs. Both 
RISS and DISS can either help or hinder pulsar detec- 
tion due to the time-variable constructive or destructive 
interference of wave fronts they cause. For a single-pass 
survey, this means that some weak pulsars may be de- 
tected because they are modulated above the detection 
threshold, and some pulsars that are otherwise bright 
enou gh may be missed beca use of scintillation modula- 
tion. ICorde£^iTaii3 (|199lD discuss in detail the detec- 
tion probability for scintillating pulsars by single-pass 
and multi-pass surveys. Since PALFA is a single-pass 
survey, it may have missed pulsars that scintillate on 
time scales comparable to or longer than the observation 
time. None of our discoveries in Table [3] appear to be af- 
fected by DISS, though RISS over long time scales may 
affect the objects with higher DMs. 
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5.2. Observation Time vs Pulsar Period 

The detectability of a pulsar as a periodic source de- 
pends on the integrated pulse flux within the observation, 
which is determined by the number of periods within the 
observation time, Np ~ Tohs/P, along with the pulse 
amplitude distribution. For the same integrated flux 
per pulse, fewer periods within a fixed Tobs mean that 
long-period pulsars are less likely to be detected than 
short-period pulsars by any method relying on detecting 
integrated flux such as an FFT-based search or a con- 
tinuum imaging survey. For FFT searches in particular, 
this effect is compounded by the fact that when A'p is 
small, the harmonics of interest are in the low-frequency 
part of the power spectrum, which is often dominated by 
non-Gaussian red noise and RFI. Hereafter we use the 
term "small- iVp bias" to refer to the combined influence 
of these effects on pulsar detectability. 

6. INTERMITTENCY MEASURE 

Both the periodicity and single pulse searches perform 
with varying efficiency depending on an object's degree 
of intermittency. In this section we present an intermit- 
tency measure method to quantify the relative perfor- 
mance of the two search algorithms, apply it to results 
from the Parkes Multibeam and PALFA surveys, and 
discuss general implications for surveys. 

6.1. Definition 

We can compare the detectability of objects by peri- 
odicity and single pulse search and attempt to narrow 
down different classes of intermittent emitters by calcu- 
lating the intermittency ratio 

r^iS/N)sp/{S/N)p^^ (8) 

for each object in our sample of PALFA and Parkes de- 
tections. iMcLaughlin fc Corded (|2003[ ) derive 

where C « 1.06 and 77 ~ 1 for a Gaussian pulse shape, 
'5'max is the maximum expected pulse intensity within 
Np = Tohs/P periods and S'^^ is a modified average pulse 
peak intensity. 

The intermittency ratio is a system-independent mea- 
sure of the efficiency of the two types of search for ob- 
jects with different degrees of intermittency. In Fig. [71 
expected values of r vs. Np are shown for an exponential 
pulse amplitude PDF and power-law distributions with 
various indices. Lower limits on r are given for inter- 
mittent sources which were not detected via periodicity 
search. In addition, we show r values for pulsars that 
were detected using both the periodicity and single-pulse 
search algorithms in the PALFA and Parkes Multibeam 
surveys. A total of 283 PALFA pulsar detections and 
255 Parkes Multibeam detections are shown in the plot, 
including multiple detections of some objects. 

6.2. Implications for Surveys 

In Fig. [71 long-period pulsars are predominantly found 
in the upper left (small Np, large r) and millisecond pul- 
sars in the lower right (large Np, small r).PSR B0525+21, 
with r ^ 7, has P — 3.7 s and DM = 51 pc cm^^. 
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Fig. 7. — The intermittency ratio r vs Np = Tohs/P for 283 
PALFA (filled hexagons) and 255 Parkes Multibeam (squares, open 
hexagons) pulsar detections. Squares denote the 10 Parkes RRATs 
with period measurements, including two which can sometimes be 
detected through their time-averaged emission. The criterion for 
inclusion of canonical pulsars was a simultaneous periodic and sin- 
gle pulse detection. PALFA points include both blind survey de- 
tections and targeted test observations of known pulsars. Arrows 
denote lower limits for PALFA and Parkes intermittent sources. 
Points with bars show the average r for a set of observations of the 
same object, with the bars denoting the maximum and minimum r 
values. Millisecond pulsar detections appear at lower right. Dotted 
lines show r for power-law pulse amplitude PDFs with indices from 
0.5 to 2.0. A ratio of 10^ was assum ed for the cutoff intensities Si 
and 52 ( McLaugh lin fc Cordes|[2003l) . The dashed line shows r for 
an exponential PDF. A single pulse search performs better than 
an FFT periodicity search for r > 1. 

which make it susceptible to both the small- iVp bias and 
diffractive scintillation and therefore it is detected more 
readily as a single pulse source. At the other extreme is 
PSR B1937-I-21, a millisecond pulsar which emits giant 
pulses but is nevertheless detected with a much higher 
S/N in the periodicity search than in the single pulse 
search (r < 0.07) because its period is small compared 
to the PALFA observation time and the normal pu lse am- 
plitude is exceptionally steady (I Jenet fc Gilll200l ). This 
indicates that the effect of Np on detectability can over- 
shadow individual properties, and in survey mode mil- 
lisecond pulsars will show up overwhelmingly as periodic 
sources. In that case, processing the data with a single 
pulse search algorithm can reveal the presence of giant 
pulses in an otherwise steady emitter. 

In-between PSR B0525-f21 and PSR B1937-I-21, from 
upper left to lower right is a relatively smooth distribu- 
tion of pulsars whose intermittency ratios are predom- 
inantly determined by the influence of Np on the de- 
tectability of all types of periodic emitters. We look 
for unusual objects among the outliers from that trend. 
There are a number of pulsars detected both in pe- 
riodicity and single pulse search in the region where 
r = 1 — 5, including PSR B0656-I-14, which occasionally 
emits bursts much brighter than the average pulse and 
would arg uably be classi fied as a RRAT if located farther 
away ( Weltevrede et aL | [200 6) . Interspersed among them 
in the region are most of the PALFA and Parkes inter- 
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mittent sources. PSR J1909+06, which was not detected 
in a periodicity search but is visible as a periodic source 
when folded with its known period, has r = 1.3, and the 
long-period PALFA intermittent sources PSR J1854-I-03 
and PSR J1946+24 have r = 1.0 and 1.8, respectively. 
Of the Parkes intermittent sources, seven have r < 3. 
PSR J1754-30 [r = 0.2) and PSR J0848-43 (r = 0.7) 
were detected through their time-averaged emission in 
low frequency follow-up observations. In the r > 5 re- 
gion we find only five objects: PSR B0525-I-21, the Crab 
pulsar, Parkes intermittent sources PSR J1317— 5759 
and PSR J1819-1458, and PALFA intermittent source 
PSR J1928-hl5. The Crab pulsar has a high intermit- 
tency ratio due to its giant pulses but is otherwise a 
relatively steady emitter as periodic emission is consis- 
tently detected even if the giant pulses are ignored. Of 
the remaining three intermittent sources none are young 
pulsars like the Crab. The pulse amplitude distributions 
of PSR J1317-5759 and PSR J 1819-1458 are describe d 
by power laws with indices ~ 1 (jMcLaughlin et 311120061 ) . 
smaller than measured giant pulse indices of ;^ 2 — 3 (e.g. 
iKinkhabwala fc Thorsettll2000l . ICordes et al1 l2004') . 

We conclude from Fig. [7] is that periodicity and single 
pulse searches should be used in combination by pulsar 
surveys since there are types of sources which are best 
detected by either method. As more intermittent ob- 
jects are found, placing them in r — iVp space may help 
us identify physically distinct populations. In addition, 
repeated sky coverage can probe intermittency at time 
scales of days to months. 

7. SINGLE PULSE DISCOVERIES 

TableOlists parameters of the seven pulsars discovered 
by PALFA via a single pulse search. In this section we 
discuss each object's properties in detail and describe the 
steps taken in addition to the PALFA processing pipeline 
in order to verify the signals' celestial origin and calculate 
the pulsar period from times of arrival of single pulses. 

7.1. PSR J 0627 +16 

PSR J0627-I-16 was discovered when nine single pulses 
were detected at a trial DM = 125 pc cm~^ (Fig. [5]). 
The period estimated from the pulse arrival times in the 
discovery observation is P = 2.180 s. What is unusual 
about this object is its very narrow pea k in S/N vs. DM 
space , indicating a narrow pulse (Corde s~fc McLaughlinI 
l2003f ). Fig. 131: shows the cumulative dynamic spectrum 
after stacking raw data chunks aligned around the five 
brightest pulses. A fit to the dispersion sweep in the 
time-frequency plane gives a best DM of 113 pc cm"'^ 
and when the raw data chunk around the brightest pulse 
is dedispersed with this DM, the FWHM pulse width is 
only 0.3 ms. 

Two follow-up observations at 0.33 and 1.4 GHz 
yielded 17 pulses in 50 minutes and 22 pulses in 72 min- 
utes, respectively. After removing bright single pulses 
from the dedispersed time series and replacing the re- 
spective time samples with a running average, periodic 
emission with P — 2.180 s was detected throughout the 
0.33 GHz observation, confirming the period estimate 
from the discovery observation at 1.4 GHz and the pres- 
ence of underlying normal emission. 

Man y of the RRATs detected by iMcLaughlin et al.l 
(|2006f ) have short duty cycles, calculated as the ratio 



between the average width of single pulses and the pe- 
riod. PSR J0627-I-16 has the smallest known duty cy- 
cle, /dc = 0.01% by this definition. However, for steady 
emitters the duty cycle is typically defined as the ratio 
of the folded pulse profile FWHM and the period. Indi- 
vidual pulses may occur in a phase window wider than 
any individual pulse and corresponding to the folded pro- 
file width. At 1.4 GHz, the individual 0.3 ms pulses of 
J0627-I-16 occur in a 8 ms window, while at 0.33 GHz 
individual pulses have widths of 0.3 — 2 ms and occur in 
a 15 ms window. The folded pulse profile at 0.33 GHz 
has FWHM of 60 ms, suggesting that bright pulses may 
be confined to a narrower window than pulses comprising 
the underlying normal emission. 

The detection of bright single pulses and a weak pe- 
riodic signal from PSR J0627+16 when data are folded 
with an accurate period estimate, along with the absence 
of a periodic detecti on in search mode weig hs in favor 
of the proposition of IWeltevrede et all ()2006l ) that some 
RRATs may have periodic emission with highly variable 
pulse amplitudes. 

7.2. PSR J0628+09 

PSR J0628+09 was discovered by detecting three 
pulses at DM = 88 pc cm^'^. The period estimated 
from the pulse times of arrival was 2.48 s. In follow-up 
observations the pulsar was also detected in periodicity 
searches as it emits on average several bursts per minute. 
The periodicity detections and the larger number of sin- 
gle pulses in subsequent observations allow ed the actual 
pulsar period of 1.241 s to be determined (jCordes et al.l 
12006) . 

7.3. PSR J 1854+03 

PSR J1854-I-03 was discovered in 2008 via a single 
pulse search performed on full-resolution survey data. 
Four pulses were detected at DM ~ 216 pc cm^"^ during 
the 268 s observation. The period P = 4.559 s was esti- 
mated by taking the smallest difference between times of 
arrival (TOAs) of two consecutive pulses and verifying 
that intervals between all four TOAs are integer mul- 
tiples of it. A confirmation observation with the more 
sensitive central ALFA beam yielded within 120 s five 
pulses whose arrival times match the estimated period. 

7.4. PSR J1909+06 

PSR J1909+06 was observed in 2006 but not identi- 
fied as a candidate until 2007 when the full-resolution 
data were searched for single pulses. Two pulses with a 
width of '--^ 1 ms and signal to noise ratio of 6 and 9 were 
detected at DM = 35 pc cm~^. The stacked dynamic 
spectrum clearly shows a dispersion sweep (Fig. |4ji). 
PSR J1909-I-06 was discovered in data from the off-axis 
beam 2 of the multi-beam ALFA receiver. Since the on- 
axis gain of the center beam is 10.4 K/Jy compare d to 
8.2 K/Jy for the other six beams (jCordes et al.ll2006l ). we 
aimed the center beam at the discovery coordinates for a 
confirmation observation. For the same integration time 
of 268 s, 8 pulses with S/N > 5 were detected in the con- 
firmation observation. We used the pulse arrival times 
to determine a period and arrived at a best estimate of 
P = 0.741 s. 
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7.5. PSR J 1919+ n 

During the 2007 discovery observation of 
PSR J1919+17, multiple pulses at DM = 148 pc cm'^ 
were detected in beam 4 of the ALFA receiver, with 
no corres ponding peri odic detection. A Fast Folding 
Analysis l(staelinlll968llKondratiev et al.ll2009D was used 
to narrow down the period to P = 2.081 s, which was 
confirmed in 2008, when the pulsar was detected as a 
normal periodic emitter with the more sensitive central 
ALFA beam. The pulsar has a double-peak profile with 
the two peaks ~ 100 ms apart, which most likely made 
it difficult to find its period from pulse arrival times 
alone. 

7.6. PSR J1928+15 

PSR J1928+15 was discovered in 2005 by detection 
of what looked like a single bright pulse at DM = 
245 pc cm~^ in a 120 s observation (Fig. [3l top). More 
detailed analysis revealed that the event was in fact 
composed of 3 separate pulses occurring at intervals of 
0.403 s, with the middle pulse being brighter by an order 
of magnitude than the other two (Fig. [3l bottom). In 
Fig. 11^ the dispersion of the brightest pulse by ionized 
interstellar gas is shown in the time-frequency plane, ev- 
idence of the non-terrestrial origin of the pulses. A fit to 
the pulse signal in the time-frequency plane resulted in a 
refined estimate of DM — 242 pc cm^'^. Despite multiple 
follow-up observations, the source has not been detected 
again. 

Given the DM of this source, 242 pc cm~^, it is unlikely 
that the non-detection can be attributed to diffractive 
scintillation. Since the three pulses are equally spaced, 
they can be interpreted as a single event seen in succes- 
sive rotations of a neutron star. This signature might be 
accounted for by an object that is dormant or not gener- 
ally beamed toward the Earth and whose magnetosphere 
is perturbed spo radically by accretion of material from 
an asteroid belt (jCordes &: Shannoiill2008D . 

7.7. PSR J 1946 +24 

PSR J1946-f 24 has DM = 96 pc cm"^ g^^^ p ^ 4 729 s 
and was discovered by detecting 4 individual pulses. The 
intervals between detected pulses in this case were all 
> 20 s. The brightest detected pulse of PSR J1946-h24 
has S/N = 29 and its dispersion sweep is visible in the 
time- frequency plane without any stacking (Fig.|4)3). 

7.8. Summary of Timing Solutions 

Single pulse arrival times collected over multiple 
follow-up observations were used to obtain partial 
timing solutions for PSR J0628+09, PSR J1909+06, 
PSR J1919-I-17, PSR J1854-t-03 and PSR J1946-I-24 
and verify their estimated periods. The periods of 
PSR J0627-I-16 and PSR J1909-h06 were verified by man- 
ually folding raw data and detecting a periodic signal. 
Finally, despite the fact that PSR J1928+15 was not 
detected after the discovery observation, the ~ 0.403 s 
intervals between the three pulses detected from that 
source differ by only ~ 2 ms, which is approximately 
half the FWHM pulse width and therefore provides a 
period estimate to that precision. We are currently tim- 
ing the six sources which were successfully redetected. 
Obtaining phase-connected timing solutions will enable 



us to compare them to other neutron star populations 
and will yield positions which will facilitate observations 
at higher wavelengths. 

8. CONSTRAINTS ON FAST TRANSIENTS 

The PALFA survey was designed to detect pulsars and 
pulsar-like objects by identifying either (or both of) their 
periodic or single-pulse emission, as discussed earlier in 
this paper. The same data may also be used to detect 
transient events from other Galactic and extragalactic 
sources. At minimum, the PALFA survey can place 
limits on the rate and amplitude distribution of radio 
transients that are shorter than about 1 second in du- 
ration. Possible source classes for short-duration tran- 
sient events, to name a few, include flare stars, mag- 
netar bursts, gamma-ray bursts (GRBs), and evaporat- 
ing black holes. GRBs have diverse time scales span- 
ning ~ 10 ms to 1000 s. Mergers of double neutron-star 
(DNS) binaries and neutron star-black hole (NS-BH) bi- 
naries are proposed sources for short bursts, and they 
may a lso emit contemporaneous radio pulses (Paczynsly 
119861 iHansen fc Lvutikovl [200l[ l. Merger rates in the 
Galaxy are estimated to be ~ 1 to 150 Myr~^ for 
DNS binaries from pulsar surveys and a factor of 20- 
30 less for N S-BH bina ries from population synthesis 
studies (E.g.. [Rantsiou et al. 2008). Event rates of de- 
tectable events depend critically on the luminosity of 
radio bursts but are likely to be comparable to the 
rate of short-duraction GRBs, and thus smaller than 
~ 1 event day~^ hemisphere"^. Another phenomenon 
that may produce rare, b right pulse s is the ann ihilation 
of prim ordial black holes ()Ree£lll977l ). iPhinnev fc Tavloil 
(|1979( l estimate the Galactic rate for such events to be 
< 2 kpc~^ y r~^. Other possibilities are discussed in 
iCordej (1200% . 

Motivated by these considerations, we discuss the con- 
straints that can be made on transient events from 
PALFA survey data obtained to date. First, we note 
that most — but not all — of the aperiodic events de- 
tected to date in PALFA data (as re ported in [7|) and 
in the Parkes Multibeam Survey by iMcLaughlin et all 
appear to be periodic when they are reobserved 
sufficiently to establish a periodicity. Additionally, all 
events detected so far in both surveys have dispersion 
measures that can be accounted for by ionized gas in 
the Galaxy (using the NE2001 model), suggesting that 
the emitting sources are closely related t o the standard 
Galac tic pulsar population. By constrast. lLorimer et al.] 
(|2007l ) recently detected a strong (30 Jy), isolated burst 
with duration ~ 5 ms that shows DM ~ 375 pc cm~^, 
too large to be accounted for by modeled foreground ma- 
terial in the Galaxy or by plasma in the Small Magellanic 
Cloud, which is somewhat near the direction the event 
was found. 

8.1. Discriminating Celestial Events from RFI 

In addition to detecting radio sources on axis in data 
streams from each of the seven beams of the ALFA sys- 
tem, it is possible to detect strong events off axis with 
low but non-zero sensitivity over almost the entire hemi- 
sphere centered on the zenith. Consequently, it is possi- 
ble to detect relatively low-rate but very bright transients 
off axis as well as the much weaker but higher rate single 
pulses that we have detected in the on-axis parts of the 
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beam pattern. 

To establish that a given event is celestial in nature, we 
need to distinguish it from terrestrial interference. For 
on-axis events, this classification process is aided by the 
directionality of the seven beams and the expectation 
that most real events are likely to be seen in three or 
fewer contiguous beams rather than in all seven beams. 
Bright events may violate this expectation, however, be- 
cause near-in side lobes are fairly large (c.f. Figure 1 of 
ICordes et al.ll2006l ). Celestial and terrestrial radio signals 
can also enter the seven beams of the ALFA receiver in- 
directly — far off the nominal pointing axis — through 
reflection and scattering off telescope support structures, 
as with any antenna. Extremely bright signals (celes- 
tial or terrestrial) may be detected after scattering into 
one or more beams, and the effective gain is approxi- 
mately that of an isotropic radiator with small gain for 
events incident far off axis from the pointing direction. 
For the Arecibo telescope, which has an intricate sup- 
port structure with many possible scattering surfaces, 
such wide-angle incidence is likely to provide multiple 
paths for radiation to enter the feed optics. Radiation 
arriving along multiple paths will constructively and de- 
structively interfere across ALFA's focal plane, so that 
the event's strength will be nonuniform in the different 
data streams. 

When radiation enters the telescope optics from wide 
angles, it is difficult to distinguish celestial from terres- 
trial signals based on the nominal pointing direction of 
the telescope and on whether they are detected in a sin- 
gle beam, a subset of beams, or all beams. Conceiv- 
ably, extremely bright events from a single radio source 
could be emitted at very rare intervals that would be de- 
tected when the telescope is pointed nominally at direc- 
tions separated by tens of degrees (but within the same 
hemisphere). Such sources and events will be exceedingly 
difficult to disentangle from terrestrial RFI. 

Of course the most powerful confirmation of a par- 
ticular event is the detection of similar events in reob- 
servations of the same sky positio n. For most of th e 
events reported in §[7] and by McLa ughlin et all (j2006l ). 
redetections have been made in this manner. They have 
also relied on establishing that their signatures in the 
frequency-time plane conform to what is expected. The 
simplest celestial signals are narrow pulses that show dif- 
ferential arrival t imes in accord with the cold-p lasma dis- 
persion law (e.g. ICordes fc McLaughlinir2003( l. Celestial 
objects may also show the effects of multipath propaga- 
tion through the interstellar medium (ISM) either in the 
frequency structure that is produced in the spectra of 
nearby pulsars or in the asymmetric broadening of the 
pulses from more distant objects. For cases where spa- 
tial information associated with the nominal telescope 
pointing direction is insufficient to constrain the celes- 
tial nature of a particular event, we must rely especially 
on a detailed study of the event's structure in the time- 
frequency plane. Indee d the time- freq u ency structure of 
the event discussed by iLorimer et al.l ()2007| ) was a key 
part of the argument for its classification as a celestial 
event. 

8.2. Simple Model for ALFA Beam Patterns 

A given reflector and feed antenna have a net antenna 
power pattern with a main lobe and near-in sidelobes 



that is similar to an Airy function with a main-lobe an- 
gular width ~ X/Da, where Da is the effective reflector 
diameter. The wide-angle part of the antenna pattern 
has far-out sidelobes that are independent of Da and 
have average gains G ~ 1, similar to that of an isotropic 
antenna. However peaks in the far-out sidelobe pattern 
can have gains significantly larger than G = 1, corre- 
sponding to special directions where scattering from the 
telescope support structure is especially efficient. For 
the Arecibo telescope, we expect sidelobes to have the 
same general properties, with the caveat that support 
structures will introduce further complexity in the far- 
out sidelobe pattern that varies as the telescope is used 
to track a source. Tracking a sky position is done by ro- 
tating the azimuth arm and moving the Gregorian dome 
along the azimuth arm, thus changing the scattering ge- 
ometry. We therefore expect the sidelobe structure to 
change significantly with azimuth and elevation of the 
targeted position on the sky. For the seven-beam ALFA 
system, the main beams by design sample different parts 
of the sky, as do the near-in sidelobes. At wider angles, 
the power pattern of each feed overlaps the others but 
with considerable variation in gain from feed to feed. 

To derive simple constraints on burst amplitudes and 
rates, we expand the gain for each feed into inner beam 
and wide-angle terms: 



G [0, (b) = G,,axP« 



+ (1-7?b): 



(10) 



where 9 and (j) are polar and azimuthal angles, respec- 
tively; Pn{0,(j)) is the power pattern for the main beam 
and near-in sidelobes; rjB is the main-beam efficiency, 
and (1 — 77b) is the average level of the far-out sidelobes. 
The boresight gain is 
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where Af. is the effective telescope area, and and 
f^MB = riB^A are the solid angles of the antenna power 
pattern and of the main beam, respectively. The system 
equivalent fiux density is 



-^sy s 2 k 



A2G 



(12) 



where k is the Boltzmann constant and Tgyg is the system 
temperature. The minimum detectable fiux density is 



sys 



^iVpoiA/I^' 



(13) 



which is a function of the event duration, VF, as well 
as system parameters. It is useful to define S'min,i as 
the minimum detectable fiux density for unit gain using 
Eq. [12] with Gmax = 1- Substituting into Eq. [13] yields, 
for two polarizations and A/ = 100 MHz, 

^..a = io"Jy (y)f^V'7S^); (14) 



w 



we note that for some directions, the system temperature 
can differ from the fiducial 30 K used. 

9. CONSTRAINTS ON EVENT RATES AND AMPLITUDES 

We now derive limits that can be placed on the rate and 
amplitudes of transient events using our simple model 
for the antenna power pattern. Let be the event rate 
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per unit solid angle from a population of sources, which 
implies a total number CTobsAil in a total observation 
time Tobs- If no events are detected, we can place an 
upper limit on the rate Cs,max for minimum flux density, 
5*111111 • We treat the response of the antenna power pat- 
tern as azimuthually symmetric about the bore axis, so 
G{9,4>) G{9), and calculate the gain as a function of 
polar angle from the bore axis. For each annulus of solid 
angle Af2, the upper bound on the rate is 



< 



1 



M Annas' 

for flux densities greater than 



5*111111 — 



Srr 



G{e) 



(15) 



(16) 



The multiplier M accounts for whether a given ALFA 
feed and receiver combination samples the same patch 
of sky or not. For the main lobes, M = 7 while for the 
far-out sidelobes, M = 1. 

We evaluate fiducial values by integrating respectively 
over the main beams and over the remainder of the power 
patterns. For a main-lobe beam width of 3.5 arc min, 
fiMB ~ lO^'^-^ deg^ and Gmax ~ 10'' °. Using the average 
gain over the main beam out to the half-power point, 
(G) «G„iax/(21n2), we have 



Sn^ i n — 



2 In 25*111111,1 
G max 



11 mjy 



10 ms 
W 



1/2 



and 



Cs, 



1 ,1,9 /461 hr 

max - ,^7^ — ~ 0.12 hr-1 deg-2 / 



Tr^MBTobs 



V Tobs 



30 K 
(17) 
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Integrating over wide angles, and using rjB ~ 0.7, we 
estimate 



Srr 



1 



VB 
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10 ms 
W 



1/2 



30 K 



(19) 

where the far-out sidelobes have an average gain of 1— jys, 
and 



1 



10 



-7.0 
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(20) 

where e « 1 is the fraction of the hemisphere actually 
covered by the far-out sidelobes. 

Results for PALFA data are shown in Fig. [51 Within 
the half-power beam width the constraints are strong on 
the flux density but weak on the event rate while the 
opposite is true for the far-out sidelobes, represented by 
the part of the curve below ( « 0.03 h~^ deg~^. The 
asympotic value of the curve corresponding to the maxi- 
mum, boresight gain differs from the fiducial value calcu- 
lated above because the latter is an average over the main 
beam. In our simple model the near-in sidelobes are not 
included explicitly so our curves overestimate 5*min at the 
corresponding polar angles. However, the limits change 
rapidly, so this does not change the salient features of the 
figure, which are the seven order-of-magnitude difference 
in constrained amplitudes and the eight order of magni- 
tudes on the event rates between the main lobs and the 
far-out sidelobes. 




10-2 ^ 





10-= 10-1 10-3 001 
if (events hr-^ deg-^) 



Fig. 8. — Constraints on the single pulse event and minimum 
detectable flux density. The curve uses a total observation time 
T = 461 h and an assumed pulse width of W = 10 ms. The 
three points, from left to right, correspond to the far-out sidelobes, 
the break between the main beam and far-out sidelobes (at ss 
0.03 h~^ deg~^), and the boresight gain. The shaded region is 
excluded by the PALFA data. 

9.1. Comparison with Other Results 

A brief summary of our results is that (a) a small 
number of events has been detected from sources that 
are consistent with their membership in the radio pulsar 
population, but with a much greater deal of modulation 
in some cases; (b) no very strong pulses have been identi- 
fied in the far-out sidelobes of the Arecibo telescope and 
(c) no bursts have been detected that are extragalactic 
in origin. 

Here we co mpar e our results with those of 
iLorimer et al.l (|2007| ). who surveyed high-Galactic- 
latitude regions that included the Magellanic Clouds 
and who discovered a single reported pulse, with 
amplitude 30 Jy. They surveyed 9 deg^ with a 
detection threshold ~ 300 mJy, about 100 times fainter 
than the detected pulse amplitude. To compare any 
two surveys, we initially assume a population of radio 
sources that is homogeneously distributed in Euclidean 
space and we ignore propagation effects that might 
limit the detectability of bursts from distant sources. 
The number of events detected TV cx fliTD^^^ for a 
survey of duration T that samples instantaneously a 
solid angle fli. For multibeam systems with iVpix pixels, 
= A^pix^^i, where fli is the solid angle of a single 
beam. We assume events have the same peak luminosity 
and pulse width W and are detectable to a distance 



The ratio of number of detected events in two 



surveys with total observation times Tobs, a and robs,b, 
and bandwidths A/q and A/h is 




(M](2l) 

^Tobs,a/ \ -^pix,a^l,a / \^bSsys,h J \A/a/ 

As pointed out bv lLorimer et al.l ()2007t ). the same sur- 
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vey that detected the 30- Jy pulse should have yielded 
many more detections if the assumptions of a volume lim- 
ited sample apply. Using Eq. [HI the lone 30 Jy pulse im- 
plies there should have been Nb/Na ~ {nia/mi,)^^'^ = 10'^ 
additional detections above a 100-timcs fainter thresh- 
old. In the 90 h of follow up observations, there should 
have been m/Na ~ [T^hs.h/Tohs.:,){ma/mbf/^ « 80 ad- 
ditional detections. Comparing the Parkes and PALFA 
surveys^, we estimate that the PALFA survey should 
have identified ~ 600 pulses above threshold. 

The PALFA survey's null result on extragalactic events 
is therefore in accord with Lorimer et al.'s non-detection 
of events weaker than the single strong event they re- 
ported. With the assumptions made, these results sug- 
gest that the non-detections of pulses weaker than the 
30- Jy pulse are inconsistent with the source population 
having a distribution that is homogeneous and isotropic. 

9.2. Caveats 

There are caveats on these results related to the as- 
sumptions about the source population and to the pos- 
sible role of propagation effects in limiting detectability. 
The minimum detectable flux density for the PALFA sur- 
vey is 30 mjy for a 5-ms pulse, implying that a pulse 
of 30 Jy strength emitted at a distance -Dcvcnt is de- 
tectable to a distance I?max given by -Dmax/-Devent ~ 
(30 Jy/0.03 Jy)i/2 ^ 32, if the lone event detected 
by Lorimer et al. is from a source at the attributed 
distance -Dcvcnt ^ 0.5 Gpc (based on assigning about 
half of the DM to the intergalactic medium), the PALFA 
survey would detect fewer pulses than estimated above 
(for constant luminosity) because the universe is not old 
enough and from redshifting of the spectrum if it is steep 
in frequency. Nonetheless, a cosmological population of 
constant-luminosity sources would extend to ^ 8 times 
the attributed distance, implying from the first form 
of Eq. [2T] that ^ 7 pulses should have been detected 
in the PALFA survey, on average, while ~ 500 weaker 
pulses should have been seen in the discovery survey 
with Parkes. The difference in yield under this alterna- 
tive scaling results from the fact that the PALFA survey 
can detect sources well beyond the cosmological popula- 
tion while the Parkes survey is comparatively shallower 
(by a factor of ten in flux density, i.e. 300 compared 
to 30 mJy), but covers a factor of 30 greater solid an- 
gle. Consideration of a broad luminosity function does 
not alter these conclusions qualitatively. However, we 
would generally expect the pulse amplitude distribution 
to extend to lower flux densities (from lower luminosities) 
given that the sole reported event was much larger than 
the threshold and that, therefore, many more pulses are 
expected compared to the constant-luminosity assump- 
tion. 

Pulse broadening from multipath propagation in ei- 
ther the ISM or the intergalactic medium (IGM) can 
also limit the numbers of events detected in a survey 
because the broadening increases with source distance. 
The maximum detectable distance is therefore dimin- 

^ For PALFA, we use 7 pixels each of size 3.5 arcmin, 461 h 
of observations, a threshold of m = 5, a system-equivalent flux 
density of 4 Jy, and 100 MHz bandwidth. For the Parkes survey 
we use 13 beams of 14 arcmin diameter, 480 h, a threshold of 600, 
Ssya = 40 Jy, and 288 MHz bandwidth. 



ished, as demonstrated for the ISM in Figure [TJ Lorimer 
et al. argued that the strong frequency dependence of 
the width of the 30- Jy pulse was consistent with multi- 
path propagation through a turbulent, ionized IGM. If 
so, sources from a greater distance would be less-easily 
detectable. In the simplest scattering geometries, scat- 
tering conserves the area of the pulse, so matched filter- 
ing yields a test statistic with S/N oc iW/rdY^'^ when 
the pulse broadenin g time Td is much larger than the in- 
trinsic pulse width (jCordes fc M cLaughlinI l2003l ). The 
30- Jy pulse therefore could have been broadened by an 
additional factor of 10'' and the S/N still would have been 
above threshold. The much longer pulses (^ lO'' x 5 ms — 
50 s) would have been difficult if not impossible to iden- 
tify in the time series owing to high-pass filtering used 
to mitigate baseline fluctuations. The high-pass filtering 
was explicit in the Parkes Multibeam data acquisition 
but is part of the post-processing in the PALFA analysis. 
Such very broadened pulses could be detected by imag- 
ing surveys with relatively short integration times like 
the NRAO VLA Sky Survey (NVSS) and future surveys 
with the Australian Square Kilometer Array Pathfinder^ 
or the Allen Telescope Array^. 

In between this extreme case of 50-s broadening and 
the reported pulse width of '-^ 5 ms is a substantial search 
volume in which more distant sources could have pro- 
duced detectable events. Broadening from sources at red- 
shifts greater than one is increased by the larger electron 
density but is lessened by the frequency redshift, so the 
net scaling of the pulse broadening with redshift depends 
on how the scattering material is distributed along the 
line of sight. 

10. CONCLUSIONS 

We have examined the current state of knowledge 
about various classes of intermittent radio-loud neutron 
stars, summarized selection effects that may affect the 
classification of intermittent sources, and presented our 
data processing methods and results from single pulse 
searches performed on PALFA survey data. One of our 
single-pulse discoveries is a relatively persistent emit- 
ter that, like many other pulsars, shows a broad am- 
plitude distribution, from which a few bright pulses 
were detected. Two sources were detected as normal 
periodic emitters when observed at a lower frequency 
(PSR J0627-I-16) and with the more sensitive central 
beam of the ALFA receiver (PSR J1919-hl7). Four ob- 
jects were most likely not detected via periodicity search 
because of their long periods (P > 2 s) compared to the 
PALFA observation time of 134—268 s, and one intermit- 
tent object (PSR J1928-f 15) was discovered by detecting 
three pulses emitted on successive rotations but it was 
not detected again despite multiple reobservations. 

Most of the PALFA sources that were discovered solely 
through the single-pulse analysis have subsequently been 
found to be sporadic but with measureable and consis- 
tent rates of detectable pulses. The Parkes survey had 
similar results. However, two sources (the Parkes source 
PSR J1839-01 and the PALFA source PSR J1928+15) 
are much more irregular even though an underlying pe- 
riodicity has been found in each case. We therefore 

^ http: //www. atnf . csiro . au/projects/askap 
^ http://ral.berkeley.edu/ata 
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distinguish between cases where the apparent sporadic 
behavior is due to a detection threshold that selects 
only the strongest puls es from a broad distribution (e.g. 
IWeltevrede et al]|2006[ ) and those where emission is truly 
intermittent. 

A major issue is related to the question, if we see tran- 
sients from a given sky direction only once, how can we 
figure out what phenomena produce them. We must 
archive radio survey data and characterize the wide vari- 
ety of detected radio transients outside the framework of 
rotating neutron stars if necessary. Archived data may 
yield new discoveries if reprocessed with improved search 
algorithms and a database of transient signals would be 
an invaluable reference as theoretical understanding of 
energetic radio events and prediction of their signatures 
improves. The PALFA survey is taking steps in both di- 
rections: raw data are stored in a tape archive and data 
processing products are indexed in a database at the Cor- 
nell Center for Advanced Computing and currently work 
is under way to make them publicly available via a web 
portal.^ 
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® http://arecibo.tc.cornell.edu 
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